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Abstract 

A classical model for the lattice distortions of Lai-xAxMnOs is derived 
and, in a mean field approximation, solved. The model is based on pre- 
vious work by Kanamori and involves localized Mn d-electrons (which 
induce tetragonal distortions of the oxygen octahedra surrounding the 
Mn) and localized holes (which induce breathing distortions). Param- 
eters are determined by fitting to the room temperature structure of 
LaMnOs. The energy gained by formation of a local lattice distortion 
is found to be large, most likely ~ 0.6 eV per site, implying a strong 
electorn-phonon coupling and supporting polaronic models of transport 
in the doped materials. The structural transition is shown to be of the 
order-disorder type; the rapid x-dependence of the transition tempera- 
ture is argued to occur because added holes produce a "random" field 
which misaligns the nearby sites. 
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LaMnOa is an insulator which undergoes a structural phase stransition at a 
Ts(x = 0) ~ 750K. The high temperature phase is believed to be cubic. The low tem- 
perature phase is approximately tetragonal, with one lattice constant rather shorter 
than the other two Several other rather small amplitude distortions also occur 
at temperatures less than or equal to Tg [0], and the structure at room temperature 
is orthorhombic. These small distortions will be ignored here. As the composition 
is varied to Lai-xAxMnOs, there are two changes. First, Ts(x) decreases rapidly 
and vanishes at x = Xs 0.2 0,11 • Second, the resistivity decreases 0. However, for 
X < Xcond ~ 0.3 and temperatures of order room temperature and higher, the material 
is still insulating in the sense that the resistivity is much higher than the Mott limit, 
and increases at T is decreased In this regime a description of the resistivity in 
terms of classical particles hopping on a lattice has been shown to be self consistent 
i- 

This paper presents a model for the x < Xcond regime and an explanation for the x- 
dependence of Tg. The physical picture is as follows: the electrically active orbitals are 
believed to be the Mn d3z2_j.2 and dx2_y2 orbitals. The mean occupancy is 1 — x [0,0. 
Because the conductivity is so low, the electrons are treated classically. It is assumed 
that a site is occupied, with probability 1 — x or empty, with probability x. The d- 
orbitals are degenerate if the local environment has cubic symmetry; the degeneracy 
is lifted by a tetragonal distortion of the local environment. Kanamori deduced 
that at X = the primary lattice distortion occurring at Tg is a staggered (tt, tt, tt) 
tetragonal distortion of the oxygen octahedra surrounding the Mn sites, driven by a 
Jahn- Teller splitting of the outer Mn d-levels; anharmonic terms in the elastic energy 
couple this to the uniform strain, producing the lattice parameter changes observed 
in early scattering experiments. Kanamori's deduction was subsequently confirmed 
by more detailed studies of the structure 0. In this paper ionic displacements will 
be explicitly included in Kanamori's model, a fit to data will be given, and the model 
will be extended to x > 0. It will be shown that the energies involved in the Jahn- 
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Teller physics are much larger than any relevant temperature, so that as long as a 
classical picture for the electrons is appropriate, a local tetragonal distortion will 
occur around each Mn site where there is an outer shell electron. At each unoccupied 
site a breathing mode distortion will occur; this will act as an effective random field 
on the staggered tetragonal distortions, and will prevent them from ordering. If the 
tetragonal distortions are not coherent throughout the lattice, they cannot couple to 
the uniform strain, and the material will remain approximately cubic. 

The model considered here is a version of the "cooperative Jahn- Teller effect", 
which has generated an enormous literature p|. Surprisingly, rather little attention 
has been paid to LaMnOs since the pioneering work of Kanamori. A Hamiltonian 
describing the orbital ordering of LaMnOs was derived from a purely electronic multi- 
band Hubbard model by Kugel and Khomskii p[ and a similar Hamiltonian has re- 



cently been derived and studied via mean field theory by Ishihara et. al. |T^, but 
atomic displacements and electron-phonon coupling have not been considered. The 
new aspects of the present paper are the explicit inclusion of the lattice degrees of free- 
dom, which allows values for the electron-phonon coupling in LaMnOs to be deduced 
from data, and the discussion of the "random field" effect of holes. 

The rest of this paper is organized as follows. In Section II the model is derived. 
In Section HI the parameters are determined by fitting structural data for LaMnOs 
to the model. In Section IV the effects of added holes are discussed. Section V is a 
conclusion. Technical details of calculations are given in several Appendices. 



II. MODEL 

In this section the energy functional is derived. The main physical assumption is 
that all degrees of freedom may be treated classically. The electrons are regarded as 
the fundamental degrees of freedom and are taken to be localized on lattice sites. In 
a classical model the hopping of electrons from site to site does not affect the energy. 
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so will be neglected. Note also that the physical mechanism primarily responsible for 
localization could be the electron-phonon coupling discussed here or the "Hubbard-U" 
effects considered by other authors ||TT|. The cause of the localization is not relevant 
to the considerations of this paper, so the electronic correlation effects need not be 
explicitly considered. 

For each fixed configuration of electrons, the phonon part of the free energy is 
minimized; the result of this minimization is the energy of that configuration of elec- 
trons. The phonons are treated in the harmonic approximation. The effect of the 
undoubtedly important anharmonic terms in the lattice energy is parametrized. Only 
some of the lattice degrees of freedom are considered. These are: 1) the vector dis- 
placements 6i of the manganese (Mn) ion on site i and 2) the scalar displacement 
u[^^ of the oxygen (O) ion along the Mn-O-Mn bond direction. Thus nf is the dis- 
placement, in the x direction, of the O atom which sits between the Mn ion on site 
i and the Mn ion on site i + x. With this restricted set of displacements one may 
discuss the Jahn- Teller distortion and the uniform strain, but not the buckling of 
the Mn-O-Mn bond or the associated rotation of the octahedra. These latter lattice 
distortions occur but, I believe, are not fundamental. 

If an electron is present on site i, it will be in a state \ipi{6) > given by a linear 
combination of the two outer d-orbitals. In the classical approximation used here the 
phase of the electron is of no significance, so one may write 

iV'i(^i) >= cos6'i|d3z2_j,2 > +sin6'i|dx2_y2 > (1) 

with < < TT. 

The lattice energy ELatt is taken to depend on the Mn-0 distance and the Mn-Mn 
distance. The unit cell i is taken to include the Mn ion at position Ri + hSi and the 
three 0-ions at positions Ri + (^|b + bu[^'*) a, where a = x, y, or z and b is the lattice 
constant. Here 6 and u are defined with reference to the ideal perovskite lattice with 
lattice constant b. In the harmonic approximation 
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(2) 

Here Ki and K2 have the dimension of energy; S and u are dimensionless. One expects 
Ki > K2. 

If an electron is present on site i, there is an electron-lattice energy given by 

EjT = A E(l - hi)[cos2^^i[vf - l(vr + vf)] + sin2^i^[vr - vf]] (3) 

i 

Here 

vf = uf - nU (4) 

and hi = if an electron is present on site i and h; = 1 if not. Finally, if there is no 
electron present on site i, all of the neighboring oxygen ions are equally attracted to 
it leading to 

Ehoic = /9A^hi[vr + vf + vf] (5) 

i 

One expects /3 ^ 1 because the force exerted on the surrounding oxygen ions by a 
Mn of the wrong charge must be much greater than the force exerted by rearranging 
the proper charge among different d-orbitals. 

For fixed values of 6i and h;, Eqs. 0,^, ^ may be minimized. The details are given 
in Appendix A. The result is most naturally expressed in terms of the parameters 



pO ^ 3 K1+K2 
^ 2K1K1+2K2 



(6) 



K - 1 



2Ki(Ki+2K2) 



as 



E = EO - hi)2 + p^hf + A cos6^i 

+ kJ: JI - hi)(l - hi+a)cOs2(^i + ^a)cos2(0i+, + (7) 
+2/?KEiahi(l - hi+a)cos2(0i+di + ^Aa) + /^^K Eia hihi+a 



Here a = ±x, y, z, iIj±z = 0, ip^x = — 7r/3, ■i/'iy = tt/S, and we have followed Kanamori 
0] by adding a phenomenological anharmonicity term with coefficient A. Cubic an- 
harmonicities exist in any realistic model of lattice dynamics. The anharmonicity is 
important for two reasons: it couples a staggered distortion to a uniform one, and 
it breaks the perfect rotational (6) symmetry found otherwise if hi = 0. The term 
added to Eq. |^ is the simplest one which accomplishes there two effects and goes into 
itself under 6' ^ + tt as required. It is derived in Appendix A. 

To each configuration of orbital occupancies {^i} corresponds an average distortion 
from the ideal cubic peroviskite structure. This may be written in terms of the oxygen 
(u) and Mn (5) displacements as 

< = Ej <(Ri - Rj)[(l - hj)cos(2ej + tPf) + /3hj] 

(8) 

= Ej - - hj)cos(2^j + ) + /5hj] 

The elastic kernels are: 

'PuK^) -Kl^k (Ki+2K 
'PsW -KlZ^k Ki+K2(l- 



(9) 
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III. FIT TO DATA 



In this section the structural information of is used to estimate model param- 
eters. The analysis is essentially that of Kanamori . A two-sublattice ordering of 
Jahn- Teller distortions parametrized by angles 6i and 62 is assumed. By fitting the 
observed atomic displacements to Eqs. and |A7|, 61, 62 and elastic constants are 



determined. By requiring that the deduced 61, 62 minimize Eq. |^ the anisotropy 
energy A is found. The experimental data for the structure are given in Appendix B 
and the mean field equations are solved in Appendix C. 

It is convenient to express the lattice distortions in terms of a staggered oxygen 



displacement Ug and a uniform strain e. By rewriting Eqs. |A2|, [A7| we obtain 



< = 5|-(cos2(^i + ij^) - cos2(^2 + ^")) 

(10) 

= (Kirk) + r)) + cos2(^2 + r))) 

In Appendix B the values e^ = -.028(-l/2, -1/2, 1) and < = 0.038(1, -1, 0) are 
derived from the data of Ref That = imphes 02 = —Oi + vr; substituting 
this into Eqs. |lO| leads to equations for 9i and A/Ki which may be solved if K2/K1 is 
given. Results are listed in Table 1. 

We now turn to the value of A. The assumption of a two sublattice distortion and 
the condition hi = implies that Eq. |^ becomes 

E = Ja[cos(60i) + cos(6^2)] + 3kcos(20i - 2^2) (H) 

By minimizing Eq. |Tl| and using ^2 = ^1 + tt we find 

A ^ -2sin(4gi) ^^2) 
K sinG^i 

Values for A/k are also listed in Table 1. 

The most important information contained in Table I is that the basic Jahn- Teller 
energy Eq is much greater than the stiffness k which orients the distortions from site 
to site. Indeed, from Eq. ^the ratio may be seen to be IK2/ (K1 + K2); as it is unlikely 
that the Mn-Mn force constant K2 > the Mn-0 force constant Ki, the ratio is less 
than 1/6. The structural transition occurring at Tg ^ BOOK in LaMnOs is therefore 
of the order-disorder type, and we may expect local distortions to persist for T > Tg. 

From Table I it is also clear that the anisotropy energy is not small, although the 
precise value depends sensitively on K2/K1. 

Now consider magnitudes of energy scales. The basic scale is Ki; this is related 
to the frequency of an oxygen bond stretching phonon Uox by 



2K1/W 

"- = V^^ (13) 



^ox 



The factor of two arises because there are two Mn-0 bonds in Eq. 0. Estimating 100 
meV> hUox ^ 30 meV and using b = 4A gives 
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300eV > Ki > 30eV (14) 

I am unaware of measurements of the phonon spectrum in LaMnOs. If, however, 
it is assumed that the phonon spectrum has a rather weak doping dependence one 
may use data optical from Lai.ssSro.isMnOa [0. The highest-lying phonon modes 
were observed at tUph ~ 70 meV. It is reasonable to assume that these are the bond- 
stretching oxygen modes of interest and that these modes are only weakly dispersive; 
thus, one may identify tUph with Uoxy and estimate Ki ^ 200 eV. 

An alternative estimate may be obtained from the mean-field approximation to 
the structural transition temperature Tg ^ 750K. This is shown in Appendix E to be 
T™^ ?a 3k, and mean-field theory overestimates Tg, so 

K > 20meV (15) 

This bound on k yields K2/Ki-dependent bounds for Ki ranging from Ki > 220 
eV (K2/K = 0.1) to Ki > 50 eV (K2/K1 = 1). Values of K2/K1 > 0.5 are most 
consistent with estimates of cUoxy < 50 meV; those of K2/K1 < 0.5 with tUoxy > 50 
meV. Combining this with the estimate Ki ^ 200 eV suggests A ~ k. This estimate 
is consistent with estimates given in a standard review [0| that typical anharmonicity 
energies are of order a few hundred kelvin. 

The estimates of Ki imply Jahn- Teller energies Eg ranging from ^ 100 meV at the 
low end (Ki ~ 30 eV) to 1 eV at the high end (Ki ~ 300 eV). The estimate Uoxy = 70 
meV implies Eq ~ 0.6 eV, slightly larger than the largest Jahn- Teller energy listed in 



a standard review |T3[. In any event, because the energy splitting between the two 
d- levels is 4Eo, it is safe to assume that at any reasonable temperature the splitting 
is frozen in. Unfortunately the splitting is difficult to measure directly because most 
methods for coupling to the d-level involve changing the valence of the Mn, which 
would bring other physics in to play. The transition should be Raman active, though. 

To summarize, it has been shown in this section that the Jahn- Teller energy of 
LaMnOs may be written 



Ex=0 = K Eia COs(2^i + 2^/'a)cOs(2^i+a + 2^a) 

(16) 

+AEiCos6^i 

If A > Tg then it is reasonable to assume that at each site 9i is near one of the three 
angles favored by the anharmonicity term, so that the system may be mapped on to a 
three-state Potts model as previously noted |§. Details are given in Appendix D. The 
result is conveniently written in a notation in which the state of site i is represented 
by a vector Q; with a 1 in one place and in the other two places; Qi = (1, 0, 0) 
implies there is a Jahn- Teller distortion with long axis along x, Qi = (0, 1, 0) means 
y and (0,0,1), z. Then 

Epotts = « E QiI"Qi+c + J' E Qi ■ Qi+2a (17) 



with a bond-direction dependent interaction given in Eq. and J' < a ferro- 
magnetic interaction between "straight-line" second neighbors which is of order k/A 
and was apparently not neglected in previous work. The second neighbor interaction 
is an approximation to the true interaction, as discussed in Appendix D. The three 
state Potts model has a transition in the x-y universality class as, therefore, does Eq. 



T6| . The second neighbor "ferromagnetic" coupling lifts the degeneracies which lead to 
subtleties in the behavior of the usual Potts model. The estimates of A suggest that 
the extreme Potts limit will not provide a good quantitative description of LaMnOs. 



IV. HOLES 

This section discusses the effects of added holes. It is clear from Eq. |^ that a 
hole on site i eliminates the Jahn- Teller distortion on site i and leads to a potential, 
/?Kcos(26'i+b + 2^/'b), which acts to orient the distortion on site i + b so that its long 
axis is along b. Thus added holes lead both to site dilution and to a field which tends 
to orient some of the neighbors of the hole in directions not compatible with long 
range order. 
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If A > (as seems to occur in LaMnOs) the angles favored by holes are compatible 
with the angles favored by anharmonicity; if A < an interesting competition arises, 
which will not be discussed here. 

In the A ^ limit the effect of added holes is particularly transparent. By 
following the derivation that led to Eq. [T^ one finds that a hole on site i produces a 
term in the energy 

ELie = /5'^ER'-Qi+b (18) 

b 

with = (—1, 1/2, 1/2) etc. Thus in this limit a hole manifestly produces a field 
which tends to orient the spins on neighboring sites. 

A Monte Carlo investigation based on Eq. |^ or on Eqs. |l^, |l^ would be desirable. 



Here simple arguments are given to estimate Ts(x). Assume the hole positions are 
uncorrelated with each other or with the configuration of Jahn- Teller orderings. To 
estimate the critical concentration, x^., of holes at which ordering vanishes, note that 
for site diluted systems (/3 = 0), Ts vanishes when the occupied sites do not percolate 
For the simple cubic lattice, the percolation threshold is about pc = 0.3 |[T^ so 
Xc(/3 = 0) = 0.7. Of course for such large values of x the model is not valid. For 
/3 — i> oo, each hole eliminates 5 sites (itself and 4 neighbors: two remain approximately 
correctly oriented), implying 1 — Sx^. = .3 or Xc(/3 oo) = 0.14. 

Alternatively, one may use mean field theory to estimate Ts(x). The fundamental 
object in mean field theory is the probability distribution P(^^) of the angle on a 
distinguished site in an effective field depending on the average values of the angles 
on the adjacent sites and on whether or not holes are present. The assumption of 
uncorrelated holes implies 

-E(e,{ha})/T 

m = E ^TTTT^^^Hl - x)^-°'> (19) 

{h,} ^li^ajj 

Here {ha} is a distribution of holes on sites adjacent to the distinguished one, Uh is 
the number of holes in that particular configuration, and 
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Z{{K})= r -e-^(^'{^'^»/^ (20) 
Jo n 

The energy may be written in terms of the average values of the cosine and sine 
on the other sublattice, c = (cos20) and s = (sin2^^) as 

E(^, {hJ)2«:Ea(l-ha)cos(2^ + 2V^,) 

(21) 

■[(c cos2^/'a — s sin2'?/'c)l — ha + Phc] 

The quantities c and s satisfy a self-consistency equation; the linearized equation 
giving Tg may be written 

f'^ d9 

c = - / — P(^) COS20 (22) 

JO TV 

The derivation and evaluation of this equation are given in Appendix E. An ana- 
lytic treatment is not simple except in the limits A ^ or A ^ oo (arbitrary j3) and 
/? — i> or /3 — i> oo (arbitrary A). For A = 0, 



^ 3k 

^ tT 



f. 41? I 



(23) 



Here the In are Bessel functions of imaginary argument i/3/Ts. 
In the A — > oo limit, Tx(/3, x) satisfies 



3k 

tT 



6 - 3e-=^^/^^- + 6e-=^^/^- , 
^ M +2e-3/3/2Ts)2 



(24) 

For (3 = 0, the x^ and higher terms vanish and Tg = 3/t(l — x) as expected for 



simple site dilution [Q. The mean field theory overestimates the x^ at which Tg 
vanishes because it does not contain the physics of percolation. As P /Tg is increased, 
the coefficient of dTg/dx increases; for /3/Tg ^ oo, Tg ^ 3k(1 — 6x), suggesting 
Xc ~ 0.16. Comparison to the percolation argument given previously suggests that 
this is an underestimate. The general result, however, of a Tg(x) which drops rapidly 
as X is increased and depends somewhat on model parameters (and so on materials), is 
in reasonable accord with data. Note however that at low T quantum effects involving 
motion of holes will become important. 
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V. CONCLUSION 



A classical model for Lai_xAxMn03 has been analyzed. It is known that doping 
on the La site changes the valence of the Mn site in such a way that the mean number 
of outer d-shell electrons on the Mn is 1-x. The holes were assumed to be classical, 
so that each Mn site is occupied, with probability 1-x, or empty, with probability x. 
The hypothesis of classical holes has been shown to be consistent with the resistivity 
at all X and T > 400K and for x < Xcond ~ 0.3 and all T ||^. 

Because the outer Mn d-orbital is twofold degenerate, a Jahn- Teller distortion of 
the surrounding oxygen octahedron which lowers the local cubic symmetry to tetrag- 
onal is may occur about each occupied Mn site, while a breathing mode distortion 
may occur around each unoccupied site. Each oxygen is shared by two Mn ions, so 
distortions on adjacent sites are coupled. The coupling was determined from a clas- 
sical harmonic approximation to the lattice dynamics. The parameters of the model 
were determined by fitting to the structural data obtained for LaMnOs. The principal 
results are: 

1) The basic energy gained in a local Jahn- Teller distortion, Eq > 0.1 eV. The 
estimate Eq ~ 0.6 eV was obtained using a phonon frequency estimated from an 
optical measurement on Lai.gsSr isMnOa. A direct measurement of the splitting 4Eo 
between the two d-levels would be desirable. The distortions are in any event well 
formed at any relevant temperature and the structural transition is to be regarded as 
an order-disorder transition, at which local Jahn- Teller distortions become spatially 
decorrelated, but do not disappear. 

2) The model describing the transition is given in Eq. ^ and may be approximated 
either by an antiferromagnetic x-y model with a modest three-fold anisotropy by or 
a three-state Potts model with an antiferromagnetic first neighbor interaction and 
a weak second neighbor interaction. Which model is more nearly correct depends 
on whether the anharmonicity parameter A is larger or smaller than the stiffness k 
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which orients the distortions. By combining an optical measurement of the highest 
phonon frequency in Lai.gsSro.isMnOs with a calculation of Tg the estimate A ~ k 
was obtained. 

3) Added holes disrupt the long range order by producing an effectively random 
field, which mis-orients nearby Jahn- Teller distortions. It would be very interesting if 
it were possible to observe directly this local misorientation. This random field effect 
was shown by various mean-field calculations to lead to a rapidly decreasing Ts(x), 
in qualitative accord with data. A Monte-Carlo investigation of the problem would 
be useful. 

The results in this paper substantiate to some degree the proposal [^,0 electron- 
lattice interaction is so strong that the high T cubic (or pseudocubic) 0.2 < x < 0.4 
phase of Lai„xAxMn03 should be modelled as a disordered array of polarons. The 
results presented here provide a basis for calculating polaron binding energies and 
mobilities, both for 0.2 < x < 0.4 and high T and for low x at all T. 
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APPENDIX A: DERIVATION OF ENERGY 

This Appendix outhnes the derivation of Eq. |^ from Eqs. 0, 0, and |^ and discusses 
anharmonic terms. Define 

ur = Ee-^''V (Ai) 

and similarly S^. Now may be decoupled from ELatt, Eq. |^, by defining 
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6^ = gives the equilibrium positions about which the Mn ions fluctuate. After 
decouphng, the relevant part of the lattice energy may be written 

ELatt = Ki^F(kJu^u% (A3) 

ka 

with 

^ i (K.+ 2K.)(1 cosu.) 

2 Ki + K2 - Kacoska 

The interaction energies are most conveniently written in terms of the variables cj^ 
defined via 

c;^ = ^e^'^-^Kl - hi)cos2(0i + tP'') (A5) 

i 

where = 0, = — vr/S and ip^ = tt/3 were introduced in Eq. 0. Combining Eqs. 
I, I, ^ gives 



EjT =AEka(l-e-'"^) 

(A6) 

Ehoic =/5AEkahk(l-e-^'^^)u^k 
The displacement u may be eliminated by writing E in terms of 

Again u = defines the average state about which the oxygen atoms fiuctuate. The 
electronic part of the energy may then be written 

E = E + /5hk)(c% + /3h_0 (A8) 

Fourier transformation yields Eq. |^ except for the term proportional to A. This term 
arises from a lattice anharmonicity of the form X^i^f. Use of Eq. |A7| yields several 
terms, of which the largest is Acos(6^). 
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APPENDIX B: ANALYSIS OF STRUCTURE 



In this Appendix the LaMnOa structural data obtained by Ellemans et. al. 
are analyzed. The magnitudes of the atomic displacements observed in LaMnOs in 
ref are somewhat greater than those reported in previous work |15|. Indeed, the 



displacements reported for LaMnOs by Ref ||T3| are very similar to those reported by 
Ellemans et. al. for Lai.gsCao.osMnOs 0. It will be assumed here that the larger 
values are most representative of the undoped material, and that the LaMnOs sample 
studied in Ref |]T3| was inadvertantly doped. 

The actual crystal structure of LaMnOs is complicated; for example every Mn-0- 
Mn bond is buckled. It is assumed here that the important quantities are the Mn- 
O bond lengths, and that the remaining distortions are subsidiary, being driven by 
rotations of the (distorted) MnOg octahedra required to fit the rigid MNOg octahedra 
into a lattice with lattice constants smaller than twice the Mn-0 distances. 

LaMnOa was found to be orthorhombic, with three unequal Mn-0 distances, which 

are 

Ux = 2.187A 

Uy = 1.905A (Bl) 
Uz = 1.956A 

Here the distances most nearly parallel to the nearest neighbor Mn- 

Mn bonds in the basic (orthorhombic a-c) plane. The x-y directions are about 45° 
rotated from the orthorhombic a-c axes, u^ is the distance most nearly parallel to the 
orthorhombic b axis. 

The z-oxygen is equidistant from the Mn above and below it; the x and y oxy- 
gen ions are not equidistant from the in-plane Mn; indeed, if one moves from one 
Mn to its in-plane nearest neighbor, the roles are reversed. We therefore assume 
that the observed Mn-0 bond lengths have been obtained from an ideal peroviskite 



15 



structure with Mn-0 distance of u*^ = (uxUyUz)^/^ = 2.0lA by composing a uni- 
form tetragonal distortion A" = — 0.112A(1, —1/2 — 1/2) and a staggered distortion 
A*^ = .15A(1,-1,0). The mean lattice constant is 4A; thus A" corresponds to a uni- 
form strain components = e^^ = .014, e™ = —.028, while the staggered distortion 
is u" = .038(1,-1,0). 



APPENDIX C: SOLUTION OF MEAN FIELD EQUATIONS 

In this Appendix details are given of the solution of the mean field equations and 
of the energetics of small deviations from the mean field solution. 
Solution: assume a two sublattice solution with 6 = 6i on one sublattice and 6 = 62 
on the other. Take A > without loss of generality and choose units in which 3k = 1. 
Write 



2ni,2 + 1 

O 



(CI) 



with 



— 7r/6 < 61^2 < 7r/6 



(C2) 



From Eq. |^ one has 

EmF = COS 

with a = A/k > 0. 
Minimizing yields 



2(ni - n2)7r 



26i — 282 



(cos65i + COS652) 

6 



2sin 



2(ni - n2)7r 



+ 25i - 2^2 



a sin6(5i 



and 



(C3) 



(C4) 



sin65i = — sin652 



(C5) 
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Equation implies either Si = —S2 or Si = S2 + (2n + l)vr/6; the latter solution 
would imply that the anharmonicity energy vanishes. Such an extremum cannot 
produce an absolute energy minimum. If S2 = —Si then Eq. may be solved. 
Define 



r ^ 1a ■ 
= Arcsm 

12 2 



'l + a2 1 



(C6) 



4a2 2a 

Then Si = —S2 = 5 if ni — n2 = 1 or —2 and Si = —S2 = —S in ni — n2 = —1 or 
2. A different formula applies if ni = n2; however, this case may be seen not to lead 
to the global energy minimum because if ni = n2 then the intersite term is positive 
unless 7t/8 < \S\ < n/Q in which case the anharmonicity term is positive. The mean- 
field energy is thus minimized by any of the six configurations with ui 7^ n2 and 
appropriate Si and S2. 

It is instructive to suppose that at all but one of the sites the angles take values 
minimizing E^p and to study the energy function Eg of the remaining angle. Assume 
the isolated site is on the "1" sublattice and ui = 0. Then 



En = 2cos 



26 - -2S 
3 . 



^cos60 (C7) 



For a < a*4/3, Eq. has only one minimum, at the 6 which satisfies the mean 
field equation. For a > a* there are three minima. For a ^ a* these occur at 



6 3a a^ 



and correspond to energies 

^ a 2n7r 2 2n7r ,^ , 

E„ = -- + 2cos^-^sm^ (C9) 

This is the expected form of the energy of a three-state Potts model with a first 
neighbor "antiferromagnetic" and second-neighbor "ferromagnetic" interactions. A 
precise mapping is discussed in Appendix D. 

Small deviations: assume that on every site the angle is close to one of the two- 
sublattice solutions; thus, if a = 1,2 
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0i = 9, + ij, (CIO) 
with 9a, given by Eqs. |Cl|, |C6| according to whether i is on sublattice 1 or 2 and 



smalL Substituting Eq. |C10| into Eq. |^ and expanding yields 

E = E^F + 5]cu2^k^_k (Cll) 
k 

The energy ul^^'^^^ depends on the quantities ni,n2 describing the possible ordered 
states. There are three independent choices ((ni, n2) = (0, 1), (1, 2), (0, 2)), each picks 
out a preferred axis a = x, y, z. We have (7k = 1 — (coskx + cosky + cosk^)) 



(0,1) 
UJC = UJ 



u = 6acos65 + 4cos ( ^ — 45) 7k 

^ ^ (C12) 



-4 



coskx — |(cosky + coskz 



Similarly a;^'^ = and w"'^ = lj^. Note that because of the relation of the angles ip 
to the physical lattice distortions, a nearly uniform variation of ip corresponds to a 
nearly staggered variation of the physical lattice distortions. For physically relevant 
values of a the gap is relatively large and the dispersion small. 



APPENDIX D: DERIVATION OF POTTS MODEL 

This Appendix gives the details of the derivation of Eq. 17 from Eq. |^ It 
is assumed A is so large that only angles near those minimizing the anharmonicity 
energy Acos66'i are allowed. Thus write 

0i = <!>? + Si (Dl) 

with (f)f one of 0x = 57r/6, 04 = 71/6 and 0z = 7r/2 and Si a small deviation. Note 
that the Jahn- Teller distortion corresponding to 0^ is Ua — |(ub + Uc). Substituting 
Eq. ^ into Eq. ^ and expanding gives 
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i = -f + I Eia COs[20r + 2^c]cos[20r+a + S^-a 

-^aEiab sm[20^ + 2^Jsin[20;^ + 2^b] 



(D2) 



cos[2V^f+, + 2V^e]cos[207+b + 2V^b] 

Note N is the number of sites in the crystal, a = A/k and (fj = 5; — 5™"^ with 
9E/95f = 0. In the large A limit the coefficient of the 6 term is large, so fluctuations 
in 6 may be neglected. 

The energy may be more conveniently written in a discrete notation. Denote the 
state on site i by the continuous variable 6i and a Q; which indicates discrete quantity 



to which of 



■'XI fyi fz 



the angle 6^ is nearest. Choose Q; to be a three component 



vector with a 1 in one place and in the other two. Q; = (1, 0, 0) means 6i is close to 
0x, Qi = (0, 1, 0) means 9i is close to 0y and Qi = (0, 0, 1) means 9i is close to (pz- The 
interaction term of order is then a 3x3 matrix I^, which depends on the direction a 
of the band connecting the two sites. One finds 



1 -1/2 -1/2 



(D3) 



-1/2 1/4 1/4 

y-1/2 1/4 1/4 y 

P, V are obtained by permuting both row and column in the obvious way. 

The term of order 1/a is a rather complicated three site interaction; however the 
important new physics of this term is the coupling it induces between sites on the 
same sublattice. To determine this coupling it is convenient to restrict attention to 
configurations (favored by the order 1 term) in which adjacent sites are in different 
Potts states, i.e. to terms in the order 1/a term in Eq. in which /5 7^ 7 or not. 
If site i is taken to be in state a = x, then the only non-constant terms are when 
a = — b 7^ X. If a = ±y then the energy is — 4/9a if both sites are in the "y" state. 
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— l/9a if both are in the "z" state, and l/9a if the two are in different states, i.e. we 
may write 



with b = ±x, y, z and 



E^'^ = -^EQi(Jb-Qi+bQi- 

ya 







-2b 



(D4) 



^b 



4-2 



^0-2 Xj 



(D5) 



etc. 



The principle effect of E'-^-' is to hft the degeneracies of the antiferromagnetic 
three state Potts model; this effect may be mimicked by a simple second neighbor 
ferromagnetic interaction with magnitude J' fixed e.g. by the requirement that it 
reproduce Eq. 



APPENDIX E: MEAN FIELD Ts WITH HOLES 

In the presence of a concentration x of holes, one expects Ts(x) = Tco(l — «x). In 
this Appendix a is derived using a mean field theory. In leading order in x one need 
only consider configurations in which one of the six neighbors of the distinguished 
site, say the one in the b direction, has a hole. From Eq. ^ one has 

Eb(^) = 2/tcos(20 + 2^b)[c cos2V^b - s sin2V^b + /?] 

(El) 

+4Ep=±i cos(2^ + 2V^b + ^)[ccos(27/^b + ^) - s sin(2^b + ^ 
Substituting Eq. |El| into Eq. Eq. expanding in c and s, rearranging, 
and discarding terms proportional to cos2^sin2'?/'b or sin26'cos2'?/;b, which will not con- 
tribute to averages of interest gives 

(E2) 

X [l - ^cos2^cos2^b - ^sin2^sin2V^t^ 
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and, using also Eq. and results of Appendix C, 



3k ri _ (-1 _ 1 1 fTT d6»„Acos6e/T-/3cos2e/T 



1 - |£^cos2^cos2V^b - ^sinesin2V^b 



xcos(26' — 2ipY^ 



(E3) 



The integrals in Eqs. p2| , [E3| may be expressed in terms of products of Bessel 
functions In(A/T)Im(/5/T). The expressions become simple when A ^ or oo or 
/3 — >• or oo, and lead after straightforward calculations to Eqs. p2|, El. 
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TABLES 



TABLE I. Values of parameters deduced by fitting mean field theory to structural data. 



K2/K1 2^}°^ A/Ki Eo/Ki X 10=^ k/Ki x 10"^ A/k 

80.9 .044 2.9 0.70 
0.1 79.2 .045 2.8 0.85 0.87 
0.3 75.7 .045 2.5 1.9 1.31 
0.5 72.3 .046 2.4 2.6 1.92 
0.75 68.2 .047 2.4 3.4 3.31 

1 64.4 .049 2.4 4.0 6.93 
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